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Abstract 

We consider the geometry of the space of Borel measures endowed with a distance that 
is defined by generalizing the dynamical formulation of the Wasserstein distance to concave, 
nonlinear mobilities. We investigate the energy landscape of internal, potential, and inter- 
action energies. For the internal energy, we give an explicit sufficient condition for geodesic 
convexity which generalizes the condition of McCann. We take an eulerian approach that 
does not require global information on the geodesies. As by-product, we obtain existence, 
stability, and contraction results for the semigroup obtained by solving the homogeneous 
Neumann boundary value problem for a nonlinear diffusion equation in a convex bounded 
domain. For the potential energy and the interaction energy, we present a non-rigorous 
argument indicating that they are not displacement semiconvex. 

Keywords: gradient flows, displacement convexity, nonlinear diffusion equations, parabolic equa- 
tions, Wasserstein distance, nonlinear mobility. 

1 Introduction 

Displacement convexity and Wasserstein distance. In |McC97] . McCann introduced the 
notion of displacement convexity for integral functionals of the form 

(//) := / U(p(x)) dx if n = p^f d , U : [0, +oo) — > E is a convex function, 
Jn 
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defined on the set ^a C (0) of the Borel probability measures in a convex open domain Q C M rf , 
which are absolutely continuous with respect to the Lebesgue measure ££ d . Displacement con- 
vexity of % means convexity along a particular class of curves, given by displacement interpo- 
lation between two given measures. These curves turned out to be the geodesies of the space 
^ c (0) endowed with the euclidean Wasserstein distance. 

We recall that the Wasserstein distance W between two Borel probability measures p,Q and [i\ 
on ft is defined by the following optimal transportation problem (Kantorovitch relaxed version) 



W 2 (no,fn) :=minW 



x-y\ dj(x,y) : 7 € r(/i ,An) 

QxSl 

where T(/io, ^1) is the set of admissible plans/couplings between /j,q and fi±, that is the set of all 
Borel probability measures on J] x O with first marginal /Uo and second marginal p,\. 
We introduce the "pressure" function P, defined by 

rr 

P(r) ■= rU'(r) - (U(r) - U(0)) = / sU"(s) ds so that P'{r) = rU"(r), P(0) = 0. (1.1) 

J 

The main result of |McC97] states that under the assumption 

P'{r)r > (1 - l/d)P(r) > 0, Vr £ (0, +00), (1.2a) 

or, equivalently, 

P(r) 



A-l/d 



is nonnegative and nondecr easing on (0, +00) 



the functional ^/ is convex along the constant speed geodesies induced by W, i.e. for every 
curve (/is) se [o,i] c ^ac(^) satisfying 

W(fJ, sl , n S2 ) = |si - s 2 | W(fJ, , m) Vsi, S2 € [0, 1], (1.3) 

the map s 1— ► ^(/i s ) is convex in [0,1]. This class of curves can be, equivalently, defined by 
displacement interpolation, using the Brenier's optimal transportation map pushing /io onto \i\ 
(see |Vil03| . for example). For power-like functions U, P 



~ xpfi ^P^ 1, p( p ) = (fL2a) is equivalent to (3 > 1 - 1/d 



U(p) = {P~ ir ^' P(p)=ff, CGID is equivalent to (3 > 1 - 1/d. (1.4) 
[plogp if 13 = 1, 

The link with a nonlinear diffusion equation. Among the various applications of this 
property, a remarkable one concerns a wide class of nonlinear diffusion equations. The seminal 
work of Otto [QtOlj contributed the key idea that a solution of the nonlinear diffusion equation 

d t p - V • (pVU'(p)) = in (0, +00) x ft, (1.5) 

with homogeneous Neumann boundary condition on d£l can be interpreted as the trajectory of 
the gradient flow of with respect to the Wasserstein distance. This means that the equation 
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is formally the gradient flow of % with respect to the local metric which for a tangent vector s 
has the form 

-V • (pVp) = s in £1 
Vp • n =0 on dQ, 



/ p\^p\ 2 dx where 
Jn 



where n is a unit normal vector to <90. Let us note here that the equation (|1.5|) corresponds via 
(TO) to 

%>-AP(p) = 0. (1.6) 

In particular, the heat equation, for P(p) = p, is the gradient flow of the logarithmic entropy 
%(p) = Jq plog pdx. Let us also note that the metric above satisfies 



inf | J p\v\ 2 dx : s + V • (pv) = in and v ■ n = on <9fi | . 



The key property of this metric is that the length of the minimal geodesic between given two 
measures is nothing but the Wasserstein distance. More precisely 

W 2 (p ,pi) = mf { f [ \v s (x)\ 2 p s (x)dxds : d s p + V • (pv) = in (0, 1) x R d , 
L JO JR d 

supp(p s ) C H, p ^ d = no, Pi^ d = Mi}- 

This dynamical formulation of the Wasserstein distance was rigorously established by Benamou 
and Brenier in [BBOOj and extended to more general situations in [AGS05j and [L07] . 

As for the classical gradient flows of convex functions in euclidean spaces, the flow associated 
with (|1.5|) is a contraction with respect to the Wasserstein distance. In |AGS05j the authors 
showed that one of the possible ways to rigorously express the link between the functional , the 
distance W, and the solution of the diffusion equation (|1.5|) is given by the evolution variational 
inequality satisfied by the measures pt = p(t, -)££ d associated with (jl.5j) : 

^W 2 (pt,v)<W(v)-W(pt) V*/e^ oc (fi). (1.7) 

A new class of "dynamical" distances. In a number of problems from mathematical 
biology [H03l IBFD061 IBD091 IDRT19] . mathematical physics [l^ll^lMlF^ISSCMlKTEnSl 
ICRS08] , studies of phase segregation |GLH IS1Q8] , and studies of thin liquid films |Ber98j , the 
mobility of "particles" depends on the density p itself. More precisely the local metric in the 
configuration space is formally given as follows: For a tangent vector s (euclidean variation) 

-V • (m(p)Vp) = s in 



(s,s) p = / m(p)|Vp| dx where 



Vp ■ n =0 on d£l 



where m : [0, +oo) — > [0, +oo) is the mobility function. The global distance generated by the 
local metric is given by 

W m,n(Po,Pi) := inf { f [ \v s (x)\ 2 m(p s (x)) dx ds : 

1 J JM. d (1.8) 

d s p + V • (m(p)v) = in (0, 1) x R d , supp(p s ) C H, P o^ d = Po, pi^ d = Pi}- 
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This distance was recently introduced and studied in [DNS09J in the case when m is concave and 
nondecreasing. Similarly to the case m(r) = r, it is easy to check formally that the trajectory 
of the gradient flow of % with respect to the modified distance W mj n solves 

dtp- V- (m(p)V [/'(/))) = in (0, +00) x ft (1.9) 

with homogeneous Neumann boundary conditions on dVt. Assuming that U"m and U"mm' are 
locally integrable, we can define in this case the function P and the auxiliary function H by 

pr pr pr 

P(r) := / U"(z)m(z)dz, H(r) := / U"(z)m(z)m'(z) dz = / P'(z)m'(z)dz, 
Jo Jo Jo 

so that 

P' = m U", H' = m' P', P(0) = P(0) = 0, 

and, at least for smooth solutions, the problem (|1.9p is equivalent to (|1.6j) . 

By means of a formal computation, detailed in Section 2, the second derivative of the internal 
energy functional % along a geodesic curve (/i s ) s e[o,i] satisfying as in (|1.3|) 



W mj n(/i sl ,At S2 ) = |si - s 2 | W m) n(/xo,/ii) Vsi,s 2 G [0,1], 
is nonnegative, i.e. > 0, if the following generalization of McCann condition (jl.2a| b) 



i. 2 

As 

holds 



P'(r)m(r) > (1 - l/d)H{r) > 0, Vr G (0, +00). (1.10a) 
It can also be expressed by requiring that 



H(r) 



m 



1-1 



l d {r) m 1 ~ 1 /rf( 



1 /' r 

■j — , / P'(s)m'{s)ds is nondecreasing in (0, +00). 



As in the case of the Wasserstein distance, in dimension d = 1 the condition (jl.lOap reduces to 
the usual convexity of U. In dimension d > 2, still considering the relevant example of powerdike 
functions U, P, m as in (|1.4p , we get 

^) = |^7 TP/3 ^f^!, m(p) = p a , P(p) = V, 7 :=« + /?-l 
[plogp if/3 = l 7 

and condition (|1.10al) is equivalent to 

6(0,1], 7>l-a/d. 

In this case the heat equation corresponds to 7 = + /? — 1 = 1 and it is therefore the gradient 
flow of the functional 

^>- (2 - a ; (1 - a) i/-^ 

with respect to the distance W mi f2 induced by the mobility function m(p) = /3 Q . 
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Another interesting example, still leading to the heat equation, is represented by the func- 
tional 



Jq y ' 
and the distance induced by m(p) = p(l — p), p £ [0,1]. Notice that in this case the positivity 
domain of the mobility m is the finite interval [0, 1], a case that has not been explicitly considered 
in [DNS09j . but that can be still covered by a careful analysis (see |LM] ). 

Geodesic convexity and contraction properties. Our aim is to prove rigorously the 
geodesic convexity of the integral functional % under conditions (jl.lOal b) and the metric char- 
acterization of the nonlinear diffusion equation (jl.9p as the gradient flow of with respect to 
the distance W mi Q (II. 8h . If one tries to follow the same strategy which has been developed in the 
more familiar Wasserstein framework, one immediately finds a serious technical difficulty, due 
to the lackness of an "explicit" representation of the geodesies for W m q. In fact, the McCann's 
proof of the displacement convexity of the functionals % is strictly related to the canonical 
representation of the Wasserstein geodesies in terms of optimal transport maps. 

Existence of a minimal geodesic connecting two measures at a finite W mj ^ distance has been 
proved by |DNS09| . However, an explicit representation is no longer available. On the other 
hand in |DS08] . following the eulerian approach introduced in |OW05j . the authors presented a 
new proof of McCann's convexity result for integral functionals defined on a compact manifold 
without the use of the representation of geodesies. Here, following the same approach of |DS08] . 
we reverse the usual strategy which derives the existence and the contraction property of the 
gradient flow of a functional from its geodesic convexity. On the contrary, we show that under 
the assumption (jl.lOap smooth solutions of (jl.9p satisfy the following Evolution Variational 
Inequality analogous to 71) 



This is sufficient to construct a nice gradient flow generated by ^ and metrically characterized 
by (|l.lip . as showed in [AGS05] and |AS06| . The remarkable fact proved by |DS08j is that 
whenever a functional % admits a flow, defined at least in a dense subset of -D(^), satisfying 
(jl.llj) . the functional itself is convex along the geodesies induced by the distance W mj ^. As a 
by-product we obtain stability, uniqueness, and regularization results for the solutions of the 
problem (jl.9p in a suitable subspace of S^i^l) metrized by W m ,n- 

Concerning the assumptions on m, its concavity is a necessary and sufficient condition to 
write the definition of W mj ^ with a jointly convex integrand [DNS09J, which is crucial in many 
properties of the distance, in particular for its lower semicontinuity with respect to the usual 
weak convergence of measures. Since m > on [0, oo) the concavity implies that the mobility 
must be nondecreasing. This is the case considered in |DNS09| . However we are also able to 
treat the case when the mobility is defined on an interval [0, M) where it is nonnegative and 
concave. It that case the configuration space is restricted to absolutely continuous measures with 
densities bounded from above by M. Such mobilities are of particular interest in applications 
as mentioned before. 




Vt G [0, +oo), Vi/ G : W m {u, p ) < +oo. (1.11) 
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Plan of the paper. In next section, we show the heuristic computations for the convexity of 
functionals with respect to W mj Q. Section 3 is devoted to introduce the notation and to review 
the needed concepts on W m) ^ from |DNS09j , Moreover, we prove a key technical regularization 
lemma: Lemma [331 Subsection 13.41 addresses the question of finiteness of W m o(/in. /ii )■ provid- 
ing new sufficient conditions on m and po, pi in order to ensure that y^ m ,Q,{po-, Pi) < +oo. After 
a brief review of some basic properties of the diffusion equation (jl.6j) . in Section 4 we try to 
get some insight on the features of the generalized McCann condition (jl.lUal b). we recall some 
basic facts on the metric characterization of contracting gradient flows and their relationships 
with geodesic convexity borrowed from |AGS05UDS08j . and we state our main results Theorems 
14.111 and 14.131 The core of our argument in smooth settings is collected in Section [SJ whereas 
the last Section concludes the proofs of the main results. At the end of the paper we collect 
some final remarks and open problems. 



2 Heuristics 



We first discuss, in a formal way, the conditions for the displacement convexity of the internal, 
the potential and the interaction energy, with respect to the geodesies corresponding to the 
distance f 1 1 . 8 [> . For simplicity, we assume that = M. d and that densities are smooth and 
decaying fast enough at infinity so that all computations are justified. 



2.1 Geodesies 

We first obtain the optimality condition for the geodesic equations in the fluid dynamical for- 
mulation of the the new distance (|1.8p . As in [B03| . we insert the nonlinear mobility continuity 
equation (11.81) 

d s p + V- (m(p)u) = in(0,l)xR d . (2.1) 

inside the minimization problem as a Lagrange multiplier. As a result, we get the unconstrained 
minimization problem 

W^(// ,/Ui) = inf sup ( / / hv s (x)\ 2 m(p s (x))dxds 
(p,v) ip ^ Jo JR d ^ 

/ [p s (x)d s ip(s,x) + m(p s (x))(v s (x) ■ Vtf)(s,x))] dxds 
+ / pi(x)ip(l, x) dx - I po{x)ip{Q,x)dx\. 

Applying a formal minimax principle and thus taking first an infimum with respect to v we obtain 
the optimality condition v = V^, and the following formal characterization of the distance 

W^Oo,Ati) = supinf { - - / / |VV'| 2 m( / 9)dxds - / / pd s ^dxds 

ib P 1 JO JR d Jo JR d 

Pi(x)i/j(1, x) dx — I po(x)ip(0, x) dx >, 

Viced J 



(3 



which provides the further optimality condition 

d s ^ + ^ m '( Ps (x))\V^\ 2 = 0. (2.2) 

We thus end up with a coupled system of differential equations in (0, 1) x M. d [DNS091 Rem. 
5.19] 

' d sP + V- (m(p)V^) =0, 
^ d s V + ^m'(p)|VVf = 0. 

2.2 Internal energy 

We use the formal equations (|2.3p for the geodesies associated to the distance (|1.8|) to compute 
the conditions under which the internal energy functional is displacement convex. If therefore 
(p s ,ip s ) is a smooth solution of (12. 3p . which decays sufficiently at infinity, we proceed as usual 
[CMV031 IVTI()3l IOW05] to obtain the following formulas: 



±&(p) = - / P(p)A^dx, 



and 

12 



^W(p) = [ (P'(p)m(p) - H(p))(M) 2 dx 
as 2 jRd 

+ / H{p){-V^j • VA^ + ^A|V^| 2 )dx 
P / ( / o)m"(p)|V / o| 2 |VV'| 2 dx. 



As usual, the Bochner formula 



-V^- VAV> + -A|Vt/>| 2 = |Hess^| 2 > 
and the fact that H(p) > 0, allow us to estimate it as 

^2&(P)> I {P'(p)m{p)-{l-l/d)H{p)){^) 2 dx- 1 - [ P'{p) m "{p)\Vp\ 2 \V^\ 2 dx. 
ds z jRd 2 y Rd 

Therefore, under conditions of concavity of the mobility m(p) and the generalized displacement 
McCann's condition (jl.lOap . the functional % is convex along the geodesies of the distance W m . 

2.3 Potential energy 

Similar heuristic formulas can be obtained for the potential and the interaction energy, as in 
|CMV03llW03] . We consider the potential energy functional 

Y(H):= [ V(x)dp, 
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with V a given smooth potential. As before, it is easy to check that the second derivative of "V 
along a geodesic satisfying (|2.3p is 



d^_ 

ds 2 



y(p) = [ m(p)m'(p) (Hess V Vi/>) ■ Vi/> dx 

JR d 

+ / m(p)m"(p)f(Vp- V^)(VF- V^) - \(Vp -VV)\V^\ 2 ) dx. 
jR d v 2 / 



This formula allows us to show that this functional cannot be convex along geodesies if m is 
not linear. Technically, the reason is the presence of the terms linearly depending on Vp. We 
present a simple example: 

Example. Let us first construct the example in one dimension. The expression for the 
second derivative of the functional above reduces to 

^(P)= I m{p) m '{p)V xx ^ x dx + \ [ m (p) m "{p) Px V x ^ x dx=:I + II 

as JR z JR 

Consider the case that V is nontrivial. Then V x ^ on some interval. For notational simplicity, 
we assume that 

V x > on [-2,2]. 

Since the mobility m we are considering is not a linear function of p there exists z > such that 
m"(z) ^ 0. Again for notational simplicity, let us assume that 



m"(z) < on 



1 3 

2' 2 



The fact that we chose V x to be positive and m" negative is irrelevant because the sign of term 
iT can be controlled by the sign of p x . Let r\ be a piecewise linear function on M: 

(I if x< -i 
r)(x) = < 1 — x if x G [—|, 5] 
I \ \tx >\. 

The fact that the function is Lipschitz, but not smooth is irrelevant; smooth approximations of 
the given 77, can also be used in the construction. Let rj e (x) = 77 (|). Let a £ Cq°(R, [0,1]), 
supported in [—1, 1], such that a = 1 on \—\, V\ and J R a(x)dx = 1. Let p e = ar] £ . Note that 
J R/ o £ dx = 1. A typical profile of p £ is given in Figure [TJ 

The test velocity (tangent vector at s = 0) we consider also needs to be localized near zero. 
A simple choice is ip £ (0) = rj e . Let p £ (s) be the corresponding geodesies given by (|2.ip and (|2.2p . 
Let us observe how, at s = 0, the terms / and II scale with e: 

I e < max m(z)m'(z) max V xx (x) —^e ~ -, 
ze[o,2] s€[-i,i] e e 

II £ < — min m(z)|m"(2;)| - min K ~ — ^. 
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Figure 1: A profile at which the potential energy is not convex. 



d 2 



Thus, for e small enough, . 
of the tangent vector gjp e (0) is 



s=0 



r V(p E (s)) < 0. Furthermore note that the square of the length 



m(p £ (0))(9^ £ ) 2 dx~- 



Thus for any AgR there exists e > such that 



ds 2 



s=0 



y{p £ {s)) + A / m(p 6 (0)) {d x ^ £ f dx < 



which implies that "V is not A-convex for any AGl. 

Let us conclude the example by remarking that it can be extended to multidimensional 
domains. In particular it suffices to extend the 1-D profile to d-D to be constant in every other 
direction and then use a cut-off. We only sketch the elements of the construction. 



We can assume that VV(0) = e^. Let p £ 



p e (xd). Let x = (aii, 



To cut-off in 



the directions perpendicular to we use the length scales 1 > i > 5 > £. Let O^s be smooth 
cut-off function equal to 1 on [— Z, I] and equal to outside of [— / — 6,1 + 5]; with | VO^s] < § and 



\D 2 9i,s\ < jp- Let pi,s,e{x) = p £ (x d )9ij{\x\). Let tpi,s,e(x) = rj e (xd)0i,s(\x\). Checking the scaling 
of appropriate terms is straightforward. 



2.4 Interaction energy 

Consider the interaction energy functional 



W{p) 



W(x - y)p{x) p(y)dxdy, 
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with W a given smooth potential. As before, it is easy to check that 

£-W(p)=[ [ m(p(x))m'(p(x))p(y) Vip(x) ■ (EessW(x -yWip(x))dxdy 
«s jR d JR d 

m(p(x)) m(p(y)) V^(y) • (Hess W(x — y)Vt^(x)) dx dy 
+ 1 I m(p(x))m"{p{x)) p{y){Vp{x) -V^(x))(VW(x -y) -V^(x))dxdy 

JR d JR d 

~\ I / m(p{x))m"(p(x))p(y)(X7p(x)-VW(x-y))\Vij(x)\ 2 dxdy. 

1 JR d JR d 

It can be demonstrated that if m is non-linear then the interaction energy is not geodesically 
convex. As for the potential energy, the reason lies in the presence of derivatives of p in the 
expression above. More precisely, in one dimension the second derivative of ^(p) reduces to 

d 2 f f 

W{p)= / m(p(x))m\p(x))p{y)^ 2 x {x)W xx {x-y)dxdy 



ds 



m(p(x)) m(p(y))il>y{y)W xx (x - y)ip x (x)) dx dy 
m(p(x)) m"(p(x)) p{y) p x (x)tp 2 x (x) W x {x - y) dx dy. 



1 

+ 2 

It turns out that the example for the lack of (semi-)convexity provided for the potential energy 
is also an example (with V replaced by W) for the interaction energy. The estimates of the 
terms are similar, so we leave the details to the reader. 



3 Notation and preliminaries 

In this section, following [DNS09] . we shall recall the main properties of the distance 
introduced in (jl.8p . For the sake of simplicity, we only consider here the case of a bounded open 
domain Q, so that it is not be restrictive to assume that all the measures (Radon, i.e. locally 
finite, in the general approach of |DNS09j) involved in the various definitions have finite total 
variation. Since we deal with arbitrary mobility functions m, these distances do not exhibit nice 
homogeneity properties as in the Wasserstein case; therefore we deal with finite Borel measures 
without assuming that their total mass is 1. 



3.1 Measures and continuity equation 

We denote by M+(R d ) (resp. M+(M d )) the space of finite positive Borel measures on M. d (resp. with 
compact support) and by M(M rf ; the space of Revalued Borel measures on R rf with finite 
total variation. By Riesz representation theorem, the space JVf^IR^; M. d ) can be identified with the 
dual space of C^(M d ;]R rf ) and it is endowed with the corresponding weak* topology. We denote 
by \v\ G M + (R d ) the total variation of the vector measure v £ M^^IR^). v admits the polar 
decomposition u = w\u\ with w S L 1 ([i/[; R d ). If B is a Borel subset of M. d (typically an open 
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or closed set) we denote by M + (B) (resp. M + (B;R d )) the subset of M + (M d ) (resp. M(M d ; R d )) 
whose measure \i are concentrated on B, i.e. /i{R d \B) = (resp. \n\(R d \B) = 0). Notice that 
if B is a compact subset of M. d then the convex set in M + (I?) of measures with a fixed total mass 
m is compact with respect to the weak* topology. If m > 0, M + (I?,m) is the convex subset of 
M + (B) whose measures have fixed total mass fJ,(B) = m. 

Let Q be a bounded open subset of W 1 . Given fio, m € M + (f2) we denote by C£^(^o — * Hi) 
the collection of time dependent measures (/i s ) s e[0,i] c M + (Q) and (f s ) se (o,i) € M(J7;]R rf ) such 
that 

1. s i— ► fj, s is weakly* continuous in M + (K d ) with h\s=o = Ho an d h\s=i = Ml! 

2. (t / s) sG (o,i) is a Borel family with Jq 1 |i/ s [(0)ds < +oo; 

3. (fx, v) is a distributional solution of 

d s [i s + V • u s = in (0, 1) x R d . 

If (/i, i/) 6 CSnl/^o — * Hi) then it is immediate to check that the total mass fi s (M d ) = /z s (f2) = m 
is a constant, independent of s. In particular, ^o(^ d ) = Hi(^- d )- 

3.2 Mobility and action functional 

We fix a right threshold M £ (0, +oo] and a concave mobility function m G C°[0, M) strictly 
positive in (0, M). We denote by m(M) the left limit of m(r) as r | M. We can also introduce 
the maximal left interval of monotonicity of m whose right extreme is 

Mf := sup jm € [0, M) : m |r 0m i is nondecreasing |. 

We distinguish two situation: 

Case A M = +oo so that m is nondecreasing and Mt = M = +oo; typically m(0) = and the 
main example is provided by m(r) = r a , a € [0, 1]. This is the case considered in |DNS09] . 
When m'(+oo) := lim r i- +00 r~ 1 m(r) = lim r | +00 m'(r) = we are in the sublinear growth 
case. A linear growth of m corresponds to m'(+oo) > 0. 

Case B M < +oo, so that < Mf < M and m is nonincreasing in the right interval [Mf, M) 
(but we also allow m to be constant or even decreasing in [0, M) with Mf = 0). Typically 
m(0) = m(M) = (in this case < M+ < M) and the main example is m(r) = r(M — r), 
or, more generally, m(r) = r a °(M — r) Ql , ao, ai € (0, 1]. 

Many properties proved in the case A can be extended to the case B, but there are important 
exceptions: we refer to |LM| for further details. Using the conventions 

a/b = if a = b = 0, , , 

a/b = +oo if a > = b, ^ ' ' 
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the corresponding action density function cf) m : R x R d —> [0, +00] is denned by 

r \w\ 2 

<t>m(P,w) = < m(p) 

[ +00 if p < or p > M. 

It is not difficult to check that, under the convention the function (f) m is (jointly) convex 

and lower semi-continuous. 

Given that m is concave and (f) m is convex, when M = +00 we can define the recession 
function ip^ : R d i-> [0, +00] (recall (|?TT]) ) 



I it? I m (V) 
(^(■U)) := lim r0 m (l, lu/r) = — - — -, m'(oo) := lim m'(r) = lim > 0. 

rt+oo m'(co) r— >+oo r— >+oo r 

We introduce now the action functional 

$ m ,n ■■ M + (R d ) x M(M d ;IR d ) -> [0,+oo], 

defined on couples of measures p G M + (IR' : '), 1/ G M(M d ; R d ). In order to define it we consider 
the usual Lebesgue decomposition p = p^£ d + p- 1 -, u = w^£ d + and distinguish the following 
cases: 

1. If the support of p or v is not contained in Q, then $ mi n(/i, v) = +00; 

2. When M < +00 (Case B), we set 

I / (j) m (p,w)dx if/i ± = 0, 1^ = 
$m,fH/^) := < Jn 

[+00 otherwise; 

notice that if <fr m ,n(p,v) < +00 then p G L°°(f2) with < p < M Jz? d -a.e. in Q and 
w G L 2 (ft;] 



3. When M = +00 and m^oo) = (Case A, sublinear growth) then 

/ </> m (p,™)dx if i/- 1 = 

$m,fH/^) : = < 

[ +00 otherwise; 

4. Finally, when M = +00 and m'(oo) > (Case A, linear growth) then we set 

( / m {p,w)dx + I ip™(w x )dp x ifv ± = w L p ± <^p^ 
[+00 otherwise. 
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3.3 The modified Wasserstein distance 

Let Q be a bounded open set. Given M^M 1 G M + (f2) we define 

W m , n (//V) :=mi^j\ m ^ s ,u s )d s y /2 : G S£ n (M° - A* 1 )} (3.2) 

= inf{Y (^m,n(/i s ,^)) 1/2 ds: G e£ n (A*° M 1 )} • (3-3) 

We refer to Thm. 5.4] for the equivalence between (|ir2J) and ([33]). YJ m ,n{fJ,°, M 1 ) = +oo 

if the set of connecting curves C£n(/z° — > /z 1 ) is empty. The following three propositions are 
proved in |DNS09| . see Theorems 5.5, 5.6, 5.7, 5.15, and Proposition 5.14. 

Proposition 3.1. The space M + (f2) endowed with the distance W mi Q is a complete pseudo- 
metric space (the distance can assume the value +oo), inducing as strong as, or stronger topology 
than the weak* one. 

Given a measure a G M + (f2), the space n [a] := {/x G M + (f2) : W m! a(/i,cr) < +00} is a 
complete metric space whose measures have the same total mass of a. 

Moreover, for every /io, Ml G M + (f2) such that W mi n(Mo, Ml) < +00 there exists a minimizing 
couple (m, v) in (I3.2p (unique, if m is strictly concave and sublinear) and the curve (Ms)se[o,il * s 
a constant speed geodesic for W m ^, thus satisfying 

Wmfiim, (jt s ) = \t- s|w m) o(/i ,Mi) Vs,t g [0, 1]. 

Proposition 3.2 (Lower semi-continuity). If Q n ,Q are bounded open sets such that £? d \ n 
weakly* converges to Sf d i^, M n G (0, +00] is a nonincreasing sequence converging to M, m n is 
a sequence of nonnegative concave functions in the intervals (0,M n ) such that 

m n /(r) > m n //(r) Vr G (0,M„») if n' < n" , lim m n (r) = m(r) VrG (0,M), 

n^oo 

and /ig, /u" are sequences of measures weakly* convergent to /io and \i\ respectively, then 

liminfW mni n>o\/^) > W m , n (Mo,Mi). (3.4) 

n— >+oo 

Proposition 3.3 (Monotonicity). Lei D f2, m > m /io, Mi ^ M + (0). T/ien i/te following 
inequality holds 

W fti sj(W,/il) < W m ^(/i ,/ii). 

Proposition 3.4. Let fc G C^°(lR d ) 6e a nonnegative convolution kernel, with f Rd k(x)dx = 1 
and supp(/c) = B\(0), and let k e (x) := e~ d k{x/e). For every /i,/io,/ii G M. + (Q) and v G 
M(H; M d ) we /lave 

$m,n £ (M*fee,"*fee) < * m) n(/i,i/) Ve>0, 
W m ,n E (Mo * fee, Mi * fee) < W m ,n(// ,/ii) Ve>0, (3.5) 

lim W m ,q £ (mo * fee, Mi * fee) = W m ,fi(Mo,Mi)> (3.6) 

e-+0 

where Q e := O + B e (0). 
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Proof. If $ m ,n(p, u) < +00 then p, v are supported in Q and |DNS09l Theorem 2.3] yields 

being p*k £ , v*k £ supported in Q £ . Notice that only the concavity of m (and not its monotonicity) 
plays a role here. A similar argument and [DNS091 Theorem 5.15] yields (|3.5|) . The limit (|3.6|) 
is an immediate consequence of (|3.4p and ()3.5p . □ 

The next technical lemma provides a crucial approximation result for curves with finite 
energy. It allows for measures to be approximated by ones with smooth, positive densities. 

Lemma 3.5. Let Q be an open bounded convex set and let {p,v) £ C£n(/io — ¥ Pi) with given 
constant mass m and finite energy & m ,n(ps, u s ) ds < +00. For every e > 0,5 £ [0,1] there 
exist a decreasing family of smooth convex sets fl £ j ft and a family of curves 
C£ne(/iQ ,<5 — > p\' S ) with the following properties 

(,f = (l-S)vi*k E + 5\*, \z : =^—J?^, p £ /(n £ )=m, (3.7) 

pf = p £ /^ d \ nE , = wf^ d \ n ., p £ ' 5 ,w £ > 5 £ C°°([0, 1] x H £ ), (3.8) 

d s p £ / + V-w £ / = ^n(0,l)xf^^ p ^>s^->o, 

\ f 1 * m ,{p(jjf a ' 5 ,^ ,<5 )ds< [\ m , n (lh,vs)ds= lim f $ m , ne (pf, uf) ds, 

c e JO JO £ > S l° Jo 

where c e := 1 + 2e. 

Proof. Let us extend (p s , u s ) outside the unit interval by setting v s = and p s = po if s < 0, 
p s = p\ if s > 1; it is immediate to check that (p, u) still satisfy the continuity equation. We then 
consider a family of smooth and convex open sets O e satisfying 0, + L>2 £ (0) C f2 e C ft + L>3 £ (0) 
and define p £ := p*k £ ,v £ := v * k £ which have smooth densities p £ s ,w £ and are concentrated 
in Q + B £ (0). We perform a further time convolution with respect to a 1-dimensional family of 
nonnegative smooth mollifiers h £ (z) := e~ 1 h(z/e) with support in [—£,£] and integral 1 

P e s ■= I p £ z h £ (s - z)dz, u £ s := / u £ z h £ (s - z) dz, 

JR JR 

with corresponding densities p £ s ,w £ . Notice that p £ _ £ = p £ '°,pl +£ = p^ and, by the convexity 
of (fi m and Jensen's inequality, we have 

<j) m (p £ s ,w £ s ) < / (f) m (p £ z ,w £ z )h £ {s - z)dz, $ mjn e(p%v e s ) < / <S> m ^{p £ s ,v £ s )h £ {s- z)dz 

JR JR 

so that, being v £ s = if s < — e or s > 1 + e, 

/ $ m ,nc{p £ s ,v £ s ) ds = [ $m,fH/^D ds < f $m,n^p £ s ,i> £ s )ds < f $ m ,n{Ps,v s )ds. 

J-e JR JR Jo 



14 



We eventually set 

A*s := /^c e s-e> v ' s := c e^c e s-e; C e := 1 + 2e 

and 

M f :=(1-<JK + 5A £ , i/^:=i£ 
It is then easy to check that all the requirements are satisfied. □ 

3.4 Couple of measures at finite W m ^ distance 

We discuss now some cases when it is possible to prove that the distance between two measures 
is finite. We already know [DNS091 Cor. 5.25] (in the case A, but the same argument can be 
easily adapted to cover the case M < +oo) that when f2 is convex and bounded 

if fii = pi££ d with ||/Oi|| L oc( M d) < M then W mj n(/io,Ml) < °c. (3.9) 



We focus on the case A, M = +oo, and exploit some ideas of [Sa08j. In order to refine the 
condition (|3.9|) . we first introduce the functions 

/ \ — 1/2 1 /'+ 00 

k m , d (r) := (V +1 / d m(r)J , K m>d (r) := - J k m , d (z)dz, r > 0. 

Observe that K m ( i is either everywhere finite or identically +oo. In particular, in the case 
m(r) = r a , K m ^ is finite if and only if a > 1 — 1/d. 

Theorem 3.6. Let Q be a bounded, open convex set of M. d . Suppose that M = +oo, m > 
0, and that K m ( i is finite. Then any two measures po,pi G M + (f2,m) have finite distance 
W mj o(/xo, pi) < +oo and the topology induced by W mj n on the space M + (S7,m) coincides with 
the usual weak* topology. In particular, the metric space (M + (f2, m), W m; n) is compact and 
separable. 

Proof. We fix an open set B with compact closure in Q, and a reference measure A = pS£ d \ B with 

A(Q) = m and < p(x) < b for Jz? rf -a.e. i in B. Since W m> Q satisfies the triangular inequality, 
the first part of the theorem follows if we show that W mj ^(A, p) < +oo for every p G M + (f2, m). 

Let r : B — > Q be the Brenier map pushing A onto p: we know that r is cyclically monotone. 
We set r s := (1 — s)i + sr with image B s C (1 — s)B + s£2 C and inverse s s = r" 1 : B s — > B, 
and v s := (r — i) o r" 1 = i — s s . It is well known that s s is a Lipschitz map with Lipschitz 
constant bounded by (1 — and that the curve p s := (r s )#A belongs to C£^(A — ► p) with 

p s = p s ££ d \ Bs , Ps = Xb s p{s s )Js, J s :=detBs s , u s = p s v s Sf d . 
Since the map r i— ► r/m(r) is nondecreasing and J s < (1 — it follows that 

<6 m ,n(M.,i/.)= / -f^|^| 2 dx= / f-ft!??^ l*-(l/)-i/l 2 p(l/)dy (3.10) 

< 4^ 'CI / Hy) - y\ 2 p(y) dy = -£f°EL w2(\, /*)■ (3.n) 



m(6(l - s)~ d ) 7 B 1 vy/ Kvy/ y m(6(l-s)~ d ) 
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Taking the square and applying (13. 3p . since 

[' ( OLziri ) 1/2 ds - ^ r (—) 1,2 z- 1 - 1 " d2 - w k m 
h U«i -«)-") l & h UwJ " Am '•' ( ' ,, 

we get the estimate 

W m>Q (A,M) < 6 1/d if m , (2 (6)W 2 (A,//). (3.12) 

A completely analogous calculation with \x := /^o (resp. /z := fi\) and /i s = ^o,s (resp. // s = ^i jS ) 
shows that 

W m ^,i- £ ,^) < feV^^fe-*) W 2 (\,m) Ve > 0, i = 0,1. 

On the other hand, taking into account that the density of /Xji_ £ is bounded by be~ d , we can 
apply (I3.12|) with /io,i- e instead of A, obtaining 

W m ,nOo,i- £ ,Ati,i- E ) < h x l d e- x K mjd (be- d )W 2 (» ,i- E ,fJ,i,i- £ ). 
Therefore, the triangular inequality yields 

W m , n (/i ,/ii) < b l ' d K m4 {be- d ) (w 2 (vo,\) + W 2 ( f i 1 , A) + e~ l W 2 ^o,i-e, 



Applying this estimate to a sequence jjL n weakly* converging to jjl (and therefore converging also 
with respect to W 2 ), since the corresponding geodesic interpolants with A [i n ,\-e converge to 
/xi_ e as n — > oo with respect to W 2 , we easily obtain 

limsupW m) n(/i n ,/i) < 2b 1 ' d K m4 (be- d )W 2 (^\)- 

n—>oo 

Since lining A m ^{be~ d ) = 0, taking e arbitrarily small, we conclude. □ 

In the next result we do not assume any particular condition on m , but we ask that fj,{ <C Jz? d 
with densities satisfying some extra integr ability assumptions. 



Theorem 3.7. Let fi be a bounded, open convex set of K and assume that M = +oo, m > 0. 

In 

Pi{x? 



If the measures fii = piJ£ d \ G M + (f2,m), i = 0, 1, satisfy 



n m(/9j(x)) 



dx<+oo i = 0,l, (3.13) 



t/ien W mi n(//o,Ati) < +oo. 



Proof. We argue as in the previous proof, keeping the same notation and observing that for 
< s < 1/2 (I3~TTD yields 

*m,n(M«, < W|(A, M). (3.14) 

When 1/2 < s < 1 we invert the role of A and \i = pJ2? d in (I3.10p obtaining 

* m ,n(/^, = / j^f^h \Hv) - y\ 2 P(y) dy (3.15) 
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where s s = (1 — s)s + si is the optimal map pushing \x onto p s and J s = detDs s 1 satisfies 
J s < s~ d . (f37T5l) then yields for 1/2 < s < 1 

<W/^ S ) < 2 d+1 / -^^(\s(y)f + \y\Ady. (3.16) 

Since the range of s(y) is //-essentially bounded, the integral in (|3.16p is finite thanks to (|3,13p . 
Integrating (I3.14|) in (0, 1/2) and (|3.16j) in (1/2, 1) we conclude that W mj n(A, p) is finite. □ 

4 Geodesic convexity of integral functionals and their gradient 
flows 

4.1 Nonlinear diffusion equations: weak and limit solutions 

We consider a 

convex density function U G W?,£(Q, M) with mU" G Li oc ([0, M)) (4.1a) 
and a pressure function P : [0, M) — ► R defined by 



P( r ) : = f m(z)U"(z)dz. 
Jo 



Let us observe that P G W,'([0,M)) is nondecreasing, continuous, and P(0) = 0. When [/ 
has a superlinear growth at +oo the corresponding internal energy functional % : D('fr) C 
M+(R d ) -» (-oo,+oo] is defined as 

<&(/*):=[ U(p(x))dx, D(<&) := {// = pJSf d G M+(M d ) : I7(p) G L 1 ^)} (4.2) 

Since J7 is bounded from below by a linear function and p has compact support, the integral 
in (14. 2p is always well defined. ^ is lower semicontinuous with respect to weak convergence in 
M+(M d ) if and only if 

U'(+oo) := lim = lim U'(r) = +oo. 

rf+oo T r|+oo 

When U'{+oa) < +oo we define the functional as 

W(ji):= [ U(p)dx + U'(+oo)p ± (M. d ), p = p^ d + p ± , 

JR d 

where p is the singular part of p in the usual Lebesgue decomposition. 

Let f2 C M d be a bounded, open, and connected set with Lipschitz boundary dQ, and ex- 
terior unit normal n. We will often suppose that f2 is convex in the sequel. We consider the 
homogeneous Neumann boundary value problem for the nonlinear diffusion equation 

d t p - AP(p) = in (0, +oo) x n, d n P{p) = on (0, +oo) x dO, (4.3) 
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with nonnegative initial condition p(0, •) = po. We also introduce the dissipation rate of & 
along the flow by 

®{p) = I |VJ y |2 dx = f cj> m {p, VP{p)) dx, VO < p G L 1 ^), P{p) G ^(fl). (4.4) 
Jn m (p) 

We collect in the following result some well established facts |Vaz07] on weak and classical 
solutions to (14.31). 



Theorem 4.1 (Very weak and classical solutions). Let us suppose that £1 is bounded and po 6 
L°°(£l). There exists a unique solution p G L°°((0,+oo) x 0) n C°([0, +oo); L 1 ^)) with P(p) G 
L°°((0, +oo) x fi) n L 2 ((0, +oo); M^ 1,2 (0)) to (|4.3p satisfying the following weak formulation 

j (pd t (-VP(p)-V(^dxdt = VC G C c °°((0,+oo) x M d ), (4.5) 

and i/ie initial condition p(0, •) = po- ^ e energy % is decreasing along the flow and satisfies 
the identity 

[ U(p(T,x))dx+ [ [ l^Ml dxdt= / U( Po (x))dx, V T > 0. (4.6) 
Jn Jo Jn m (PJ Jn 

T/ie map po i— » S^po := /°(^> •) cara ^ e extended to a C° contraction semigroup S = S(P, O) in the 
positive cone o/L 1 (J7), whose curves Stpo are also called "limit L 1 -solutions" of (|4.3p . and ii 
satisfies 

ess info po < S t po < esssup Q p . 

If moreover U, m G C°°(0,M) ; [/ is uniformly convex, 0, is smooth and inf q po > 0, t/ien 
p G C°°((0, +oo) x O) and is a classical solution to (|4.3p . 



Let us briefly discuss here two useful lemma, whose proof follows from a standard variational 
argument. 

Lemma 4.2. If po, U(po) G L 1 (17) £/ien the limit L l -solution p = S(po) satisfies P(p) G 
£f oc ([0, +oo); W ' (^)), fie weak formulation (|4.5p . and fie energy inequality 

[ U(p(T,x))dx+ ( T [ |VP , ( ^ |2 dsd*< / U(po(aO)dx. (4.7) 
Jn Jo Jn m (p) Jn 

Proof. Let us first show that we can find a constant C depending only on P, u := J£ d (Q), 
m = f n pdx, and the constant c p in the Poincare inequality for such that 

\\P(p)\\ L i(n)<c(l + \\VP(p)\\ L i (Q) ) VpeL\n), ^Pdx = m, P(p) sW 1 ' 1 ^). (4.8) 

In fact, setting p := f n P(p) dx and £ := ^f d ({x G : -P(p) > p/2}) Poincare and Chebyshev 
inequality yield 

^p(co -£)< j \P( P ) - p| dx < c p f \VP{p)\ dx, l -p < P(m/£), 
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so that if £ > ui/2 we get p < 2P(2m/cj), whereas if £ < uj/2 we obtain p < Alu~ 1 c p f n |VP(p)| dx. 

If now p t = Stpo is the L 1 (r2)-limit of a sequence p n ,t = Stp n ,o of bounded solutions with 
U(p n fi) — > U(po) in L x (r2) as n | +00, from the uniform bound f)4.6fl we obtain for every 
bounded Borel set 7 C (0, +00), every P C fi, and every nonnegative constants a, b such that 
m(r) < a + br, 

JJ B \VPM\6x*< l|m(p„)||^ xB) ( fj^f dxdt) V2 
< C||a + 6p n ||^ 2 (TxS) . 

Taking T = (0,T),B = and applying (|4,8|) . we obtain a uniform bound of the sequence 
P(p n ) in L^C^T; W 1 ' 1 ^)); since p n converges to p in L^^T) x ft), we obtain that VP(p n ) 
is uniformly integrable and therefore it converges weakly to VP(p) in L 1 ^ ^) x 0). It follows 
that P{p) € L 1 (0, T; VI /1,1 (0)) and we can then pass to the limit in the weak formulation (I4.5D 
written for p n , obtaining the same identity for p. The inequality (|4.7p eventually follows by the 
same limit procedure, recalling that the dissipation functional (|4.4|) is lower semicontinuous with 
respect to weak convergence in L 1 ^). □ 

The following stability result is used in the sequel; its proof is an easy adaption of [Vaz07[ 
Prop. 6.10]. 

Proposition 4.3. Let f] n C M. d be a decreasing sequence of open, bounded, convex sets converg- 
ing to f2 and let S n = S(P,£l n ), S(P,Q) be the associated semigroups provided by Theorem \4-l\ 
If (after a trivial extension to outside Q n ) p$ € L 1 (O n ) is converging strongly in to 
po £ then S™(Po) — > St(po) in the same L 1 sense, as n j +00 for every t > 0. 

4.2 The generalized McCann condition 

We assume that P'm' g L\ oc {[Q,M)) and we introduce a primitive function H of h := P'm' = 
[/"mm', 

H(r) := H + / P'(z)m'(z) dz for some F > 0. (4.9) 
Jo 

When the dimension d is greater than 1, we assume that inf re (o,Af) -^( r ) = 0> this means that the 
locally integrability assumption on h cannot be avoided, as well as, in the case when M < +00, 
its integrability in (0, M). In the case M = +00 we can simply choose Hq = in (|4.9p . In the 
case M < +00 we can choose 

P := (J P'm'dx)" < +00. (4.10) 

Remark 4.4. Notice that in the most common case when m'(0+) = lim r |o r~ 1 m(r) > 0, the 
local integrability of h in a right neighborhood of implies the local integrability of m U", as 
we already required in (]4.1ap . and the lower boundedness of P. When the space dimension is 
1 all these restriction can be removed: we comment on this issue in the next remark 14.161 If 
M < +00, P is locally Lipschitz near and M, and m(0) = m(M) = 0, then we get Po = so 
that P and m should satisfy the compatibility condition J P'm' dr = 0. 
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Definition 4.5. Let U, P, H and m be defined in the interval (0, M) according to (|4.1al b) and 
(|4.9p . We say that the energy density U and the corresponding pressure function P satisfy the 
d- dimensional generalized McCann condition for the mobility m, denoted by GMC(m,d), if for 
a suitable choice of Hq 

U"(r)m 2 (r) = P'(r)m(r)>(l-l/d)H(r)>0, Vr€(0,M), (4.11a) 

or, equivalently, 

H(r) 



m 



l-i 



/d( r ) m l-i/d( 



i r 

- — TTj— - / P'(s)m'(s)ds is nondecreasing in (0, +oo). 
l ~ l l d (r) Jo 



We collect in the following remarks some simple properties related to this definition. 
Remark 4.6 (Elementary properties). 

1. (Linear mobility) (|4.11ajl is consistent with the usual McCann condition (jl.2aj> in the linear 
case of m(r) = r. 

2. (Dimension d = 1) As in the case of McCann condition, in space dimension d = 1 (|4,lla|) 
is equivalent to the convexity of U or to the monotonicity of P. 

3. (Local boundedness of U when d > 1). In dimension d > 1 the energy density function 
U is bounded in a right neighborhood of (and in a left neighborhood of M, in the case 
M < +oo). Since U" = P'/m the property is immediate if m(0) > 0. If m(0) = then 
m'(0) > and therefore P is bounded around and the formula 

U(r) = U(r Q ) + U'(r Q )(r - r ) + f ° (z ~ r) P\z) dz, r E (0, r ], 

Jo m \ z ) 

shows that lim r j^o U(r) < +oo. 

4. (Constant mobility) When m(r) = c > (|4.11ap is still equivalent to the convexity of U. 

5. (The power-like case) In the case of P(r) = r 7 (7 = a + (3 — 1 if U(r) = r@) and m(r) = r a , 
(14.1 lap is satisfied if and only if 

a 

7 >1--. (4.12) 

6. (The case P(r) = r) It is immediate to check that the couple (r, m) always satisfies (|4.11ap : 
it corresponds to the entropy function U m whose second derivative is m _1 . After fixing some 
ro € (0, M) (the choice ro = is admissible if m^ 1 is integrable in a right neighborhood of 
0), we obtain 

f r r — z 

U m {r) := / — — dz, P m {r) = r - r . 
Jr m(2) 

7. (The case of the logarithmic entropy) U(r) = rlogr satisfies GMC(r a ,d) if and only if 
7 = a > d/(d + 1). 
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8. (Linearity) If Pi and P2 satisfy GMC(m,d) then also a\P\ + a^Pi satisfies GMC(m,d), 
for every a\,a.2 > 0. Analogously, if P satisfies GMC{m\,d) and GMC(tr\2, d) then P 
satisfies GMC{ct\m\ + 021112, d). In particular, if P satisfies GMC(m,d) then P(r) + ar 
satisfies GMC(m + (3, d) for every a, (3 > 0. 

9. (Shift) If M = +00 and P satisfies GMC(m,d) then P satisfies GMC(m(- + a),d) and 
P(- — a) satisfies GMC(m, d), for every a > 0. 

The next two properties are more technical and require a detailed proof. 

Lemma 4.7 (Smoothing). Lei us assume that P satisfies GMC(m, d) and let us fix two constants 
< M 1 < M" < M. Then there exists a family P^rn^, i] > 0, with smooth restriction to 
[M',M"] such that P^ > P is strictly increasing, > m is concave, P v satisfies GMC(m n ,d) 
(in [0,M"]), and P^m^ converge uniformly to P, m in [M',M"] as rj [ 0. Moreover, if P' is 
locally integrable in a right neighborhood o/O, then we can choose M' = 0. 

Proof. When M' > it is not restrictive (up to choosing a smaller M') to assume that M' is a 
Lebesgue point of the derivative of P. Let H be as in (|4.9p and let us set rh^(r) := m(r) + 77, 
P v (r) = P(r) + rjr, 

H v (r) = H + 77m (0) + rf + / P'(r)m' (r) dr = iT(r) + r/m(r) + r/ 2 > r/ 2 > 0. 

•/ 

By the previous Remark (points 6 and 8) P v satisfies GMC(m,, d) and moreover 

P^m v - (1 - l/d)fl, = P'm - (1 - l/d)# + ^(r/ + m) > % m + v ) > ^/^ ( 413 ) 
By choosing a family of mollifiers h$, S > 0, with support in [0, 5], we introduce the functions 



P v ,s(r) 



P V (M') + $ T M , *h s ds if r > M' . 

P v (r) if r<M', 

' m v (M') + f M , m' n *h 5 ds if r > M', 

m^(r) if r<M', 



which are smooth in [M',M"] and satisfy the requested monotonicity/concavity conditions. 
Since P^ § converges to P; ; in L} oc (Q,M'] and m„ ,5 is uniformly bounded and converges pointwise 
a.e. to as 8 — > 0, we conclude that the corresponding continuous functions converge 
uniformly to H v as 5 [ 0. By (I4.13|) . we can find a sufficiently small 5 = 5 V depending on r\ such 
that 

P^ 5v rn' vA >(l-l/d)H v ,s v >0. 
A standard diagonal argument concludes the proof. □ 
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Lemma 4.8 (Minimal asymptotic behaviour). When d > 1 and M = +00, the function 
Pmm( r ) '■= Jo m(z) _1 / rf dz satisfies GMC(m,d) and provides an (asymptotic) lower bound for 
every any other P, since for every ro > there exists a constant cq > such that 

P'(r) > c P^ in (r) = com(r)' 1/d , U"(r) > com^)" 1 " 1 ^ fora.e.r>r . (4.14) 

Proof. In fact f(r) := P'm satisfies 

f'T 

f(r) > (1 - 1/d) (F(r ) + / /(r)m'(r)/m(r) dr) 

Gronwall Lemma then yields f)4. 14|) with cq := (1 — l/<i)//(ro)m(ro) 1 / rf_1 . □ 

Notice that in the case m(r) = r a we obtain the functions P m i n (r) = cr 10 with exponent 
70 = 1 — a/d, which is consistent with (|4.12|) . The corresponding energy density functions are 
then C/ m i n (r) = cr 2 - a ( 1 + 1 / rf ) : [ n particular, when a < d/(d+ 1), all the energy functions have a 
super linear growth as r f 00. 

Remark 4.9 (A sufficient condition). It is possible to give a simpler sufficient condition than 
(I4.11ap . at least when mil" is integrable in a right neighborhood of and M = +00: if 

the map r 1— > m (r) is positive and nondecreasing in (0,+oo) (4.15) 

then (I4.11ap is satisfied. In fact, assuming U smooth for simplicity, (|4.15p is equivalent to 

< mWp" - (1 - l/c^m^-WP'. 

Multiplying the inequality by m 1 ~ 1 / d and integrating in time we get (|4. 1 lap . Condition (|4,15p 
gives the same sharp bound (|4.12|) in the power case. 

4.3 The metric approach to gradient flows 

We recall here some basic facts about the metric notion of gradient flows, referring to [AGS05] for 
further details. Let (D,W) be a metric space, not assumed to be complete, and let "V : D{Y) — > 
(—00, +00] be a lower semicontinuous functional. A family of continuous maps St : D — > D, 
t > 0, is a C°- (metric) contraction gradient flow of "V with respect to W if 

S t +h(u) = S h (S t (u)), limS t (u) = S (u) = u Vu G D, t,h>0, 
^W 2 (S t (u),v) - ^W 2 {u,v) < t(y(v) - y(S t (u)f) Vt > 0, u e D, v e D(y). (4.16a) 
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Thanks to [DS081 Prop. 3.1], conditions (jlTTBk b) imply 

St(D) C D[Y) Vi > and the map t i— > ^ / (St(u)) is not increasing in (0, +co), 

-w 2 {s t (u),v) + r(s t (u)) < vti e a u e z>(r), t > o, (4.17) 



2 dt 



W 2 (S 4l (u),S i0 (n)) <2(ti-to)(r(Si u)-ri„f) VuGD(n 0<i <*i, (4.18) 
W(S t (u), S t («)) < W(u, u) V u,v £ D, t > 0. 
In (|4.17p we used the usual notation 

— C(t) =hmsup . 

dt h ^ + /i 

for every real function £ : [0, +00) — > R. 

The following approximated convexity estimate |DS08l Theorem 3.2] plays an important role in 
the sequel. 

Theorem 4.10 (Approximated convexity). Let us suppose that S is metric contraction gradient 
flow of "f with respect to W according to (|4.16a .b) and let s 1— » u s G D , s G [0, 1], 6e a Lipschitz 
("almost" geodesic) curve such that u$,u\ € D{Y) and 

W(u r ,u s ) < L\r - s\ Vr.se [0,1], L 2 < W 2 (« , «i) + <5 2 - (4.19) 

27ien for every s G [0, 1] and t > 0, we /iawe 

^(s t («.)) < (1 - s)r(«o) + s y( Ul ) + ^^ 2 . 

In particular, if u s is a minimal geodesic, i.e. (|4.19p /ioWs iwitt 5 = 0, then 

y(u s ) < (1 - a)r(«o) + sY( Ul ) Vs G [0, 1]. 

4.4 Main results 

We state our main result about the generation of a contractive gradient flow of % with respect 

to W m> Q. 

Theorem 4.11 (Contractive gradient flow). Let us assume that £1 is a bounded, convex open 
set, and the functions U,P,H satisfy the generalized McCann condition GMC{m,d). For every 
reference measure a G M + (f2) with finite energy &(cr) < +00 the functional % generates a 
unique metric contraction gradient flow S = , m,£l) in the space 

D:={p£ M + (fi) : [i < J2? d | n , W m ,nOu, a) < +00, <gr(/i) < +00} 

endowed with the distance W m) Q. Moreover 8 is characterized by the formula §tfio = Pt^ d \^ 
where pt = Stpo is a limit L 1 -solution of (|4.3p . 
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When m satisfies the finiteness condition of Theorem 13.61 (in particular m(r) = r a with 
a > 1 — 1/d) we obtain a much more refined result, which in particular shows the continuous 
dependence of S on the weak* topology. 

Corollary 4.12. Under the same assumptions on £1, U,P of the previous theorem, if moreover 
M = +oo and m satisfies the finiteness condition of theorem \3.6\ then the semigroup 8 can be 
uniquely extended to a contraction semigroup on every convex set M + (0, m), which is continuous 
with respect to the weak* convergence of the initial data. If U has a superlinear growth, then 
§t(A*o) = <C -^Iq f or every t > and pt is a weak solution of (14, 3p according to (|4.5|) , 

We conclude this section with our main convexity result. 

Theorem 4.13 (Convexity). Let us assume that Q is a bounded convex open set, and the 
functions U,P,H satisfy the generalized McCann condition GMC{m,d). For every po, pi € 
M + (0)nD(^') with finite distance W mj o(/xo, fix) < +oo there exists a constant speed minimizing 
geodesic for W mj n, /U : [0, 1] — > M + (0) connecting /j,q to fix such that 

W(ij s )<sW(»x) + {1-s)W(h ), Vse[0,l]. (4.20) 

Remark 4.14 (Weak and strong convexity). When (|4.20p holds for all the (constant speed, 
minimizing) geodesies, the functional % is called strongly geodesically convex. When m is 
strictly concave and has a sublinear growth (or M < +oo) then every two measures with finite 
W mj ^-distance can be connected by a unique geodesic |DNS09l Theorem 5.11], so that there is 
no difference between strong or weak convexity and (|4.20|) yields that the map s i— > fy{p s ) is 
convex in [0, 1]. 

Remark 4.15 (Absolutely continuous measures). Even when geodesies are not unique, the 
proof of Theorem 14.131 shows in fact that (|4.2U|) is satisfied by any geodesic p s with fi s <C Jz? d 
for every s € M. d , which surely exist if U has a superlinear growth. Along this class of geodesies 
we still obtain that the map s *— > ^(/i s ) is convex in [0, 1]. 

Remark 4.16 (The one-dimensional case). When the space dimension d = 1, then the general- 
ized McCann condition GMCim, 1) reduces to the usual convexity of U. In this case, a simple 
approximation argument shows that we can cover also the case of functions U which are not 
bounded in a right neighborhood of (and in a left neighborhood of M, if M < +oo) and the 
integrability assumptions on U"m of (|4,lap and on [/"mm' of (|4.10p can be dropped. 

5 Action inequalities in the smooth case 

In this section we assume that Vt is a smooth and bounded open set. We consider a smooth curve 

Us := Ps^ d \ w peC°°([0,l]xO), Q<m <p<mx<M, /x fl (fi) = m, s € [0, 1]. (5.1) 

We also assume that P and m are of class C°° in [mo, mi]. We consider the semigroup S = 
S(P,Q) defined by Theorem 14.11 and we set 

p s y.= Ps ±££ d \ p a , t (.) =p(s,t,-) :=S stPs , »€[0,1], t > 0. (5.2) 
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Classical theory of quasilinear parabolic equation shows that p € C°°([0, 1] x [0, oo) x Q) n 
C°°([0,1] x (0,+oo) x O). 

Since the semigroup St preserve the lower and upper bounds on p and j n d s pdx = 0, for 
every (s,t) € [0, 1] x [0, +oo) we can introduce the unique solution £ S)t = C( s >^ ') £ C°°(Cl), of 
the uniformly elliptic Neumann boundary value problem 

-V • (m(p)VC) = d s p in O, 

VC • n = on dQ, (5.3) 
/ n C(x)dx = 0. 

It is easy to check that £ depends smoothly on s and t. Notice that (|5.3[) is equivalent to 

f m(p)VC- Vr/dx= / S^da; Vr? e C 1 ^). (5.4) 

By construction, for every t > the curve s i— ► (/i S) t, v s ,t) with i/ a ^ := m(/? Sjt )VC Sj j^f rf | fi belongs 
to C£.(po — * /^lt) an d its energy can be evaluated by integrating the action 



A s ,t ■= ^m,n(Ms,f = / m (/°s,i)|VC s ,t| 2 da: 



with respect to s in the interval [0,1]. The integral provides an upper bound of the W mj Q- 
distance between po and p,± t t = pi fj£? , which corresponds to the solution of the nonlinear 
diffusion equation (|4,3p with initial datum px- As it was shown in [DS08j . evaluating the time 
derivative of the action A S) t is a crucial step to prove that (|4.3p satisfies the EVI formulation 
(|4.16ap . Next lemma, which does not require any convexity assumption on collects the main 
calculations. 

Lemma 5.1. Let p s ,p s j, and £s,t be as in (|5.ip . (|5.2|) . and (|5.3p . T/ien /or even/ (s,t) E 
[0, 1] x (0, +oo) we have 

2Ql A s,t = ^2 J^WCs,t\ 2 m{p s ,t)dx = - J VP(p., t ) • VC s , t d:r 

- ^ (VGomM/V) - ^(p s ,t))(AC s ,t) 2 + H(p s ,t)\D 2 t s , t \ 2 ) dx (5.5) 

+ s / P'(p s> t)m"(p s ,t)\Vp s>t \ 2 \VCs,t\ 2 dx + ls [ H(p s ^V\VC s ,t\ 2 -ndJ^ a 
Jn 1 Jan 



od-l 



where H is defined in (14, 9p . 

Proof. For keep the notation simple, we omit the explicit dependence of p, £ on s,t. By the 
definition of p we easily get 

d t p = sAP(p) and s VP{p) -Vi]dx = - d t pr]dx VijeC 1 ©. (5.6) 

Jn Jn 
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Further differentiation with respect to s yields 

- / VP(p) • Vrjdx + s / d s P(p)Arjdx = / d st pr]dx (5.7) 
,/n Jn Jn 

for all G C 2 (f2) with V77 • n = on <9il. On the other hand, differentiating (|5.4p with respect 
to i we obtain 

/ m(p)<9 t VC • Vr?dx = / d st pr]dx- / <9 t m(p)VC ■ V77 dx Vt/ G C 1 ^). (5.8) 
Jn Jn Jn 



>n Jn Jn 

The time derivative of the action functional is 
d_l 
dtl 
EH) f 



(57) 



VC| 2 m(p)dx = ^ / t m(p)|VC| 2 dx + / d t VC • VCm(p) dx 

/ d st pQdx- \ [ ^m(p)|VCj 2 dx 
Jn z Jn 

[ VP(p)-VCdx + s f d s P(p)A(dx-}- [ d t m(p)\V(\ 2 dx = I + 11 + III. (5.9) 
Jn Jn 2 7n 



We evaluate separately the various contributions: concerning the second integral II we introduce 
the auxiliary function G 



and we get 



G(r) := P'(r)m(r) - H(r), so that G'(r) = P"(r)m(r), (5.10) 

II = s f d s P(p)A(dx = s [ P'(p)d sP A(dx 
Jn Jn 

^ s [ P'(p)m(p)VA( ■ VCdx + s f A( P"(p)m(p)Vp ■ V( dx 
Jn Jn 

™s f P'(p)m(p)VAC • VCdx + s / AC VG(p) • VC dx 
in Jn 

= a / H(p)VA(- VCdx-s [ G(p)(A() 2 dx. 
Jn Jn 



The third integral of (|5.9f) is 

/// = _! /" m >(p)d tP \v(\ 2 dx = 

1 Jn 

[ P'(p)m // (p)|V /0 | 2 |VC| 2 dx + ^ / P>(p) m >(p)V\V(\ 2 -Vpdx (5.11) 
io z Jn 

[ P'(p) m "(p)\Vp\ 2 \VC\ 2 dx+ S - [ V|VC| 2 -ViJ(p)dx. 
in z Jn 

A further integration by parts in the last integral and the Bochner formula 

VC • VAC - ^A|VC| 2 = -\D 2 (\ 2 
yield (|51)]) . □ 
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Corollary 5.2. Under the same notation and assumptions of Lemma \5.1\ if£l is convex and U 
satisfies the generalized McCann condition GMC(m,d) < \4-l b), then 



\wt A °* = m i I^mIMpm) dx < = ~ J a u(p. tt ) dx, 



(5.12) 



and 
1, 



W^o(/i M ,/i )+t^(/il > t) < 



\f A s>t ds+ [ W(ji ltT )dT< \ f A afi d8 + W(jio). (5.13) 
2 Jo Jo 1 Jo 

Proof. We determine the sign of the terms in the right hand side of (|5.5p thanks to (|4.11a .b) 
and the convexity of 0. Recalling that |-D 2 C| > h(^C) 2 an d H > we obtain that the second 
integral in the right-hand side of (15. 5p is nonpositive 

- J ((P'(p)m(p) - H(p))(AC) 2 + H(p)\D 2 (\ 2 ) dx ™ 0. (5.14) 

Since P is increasing and m is concave, we have P'(p)m"(p) < which yields 

f P'(p)m"(p)|Vp| 2 |VC| 2 dx < 0. (5.15) 
Jn 

Since H is nonnegative, the smoothness and convexity of O and the smoothness of £ yields, see 
[Gri85llOtnillGSTn8j . 

V|VC| 2 -n<0, on dCl, [ H(p)X7\V(\ 2 ■nd-J^ d ~ 1 < 0. 

Jan 

Combining (I5.14|) . (I5.15|) and (I5.16|) . (15. 5|) yields the inequality 

/ |VC| 2 m(p)dx < - f V(-VP(p)dx. 
Jn Jn 



a_i 

dt2 



(5.16) 



(5.17) 



On the other hand 




ds 



[ U(p)dx= [ U'(p)d s pdx^ [ m(p)VU'(p) ■ V(dx = [ X7P(p)-V(dx, (5.18) 
Jn Jn Jn Jn 



so that (15.120 follows by (|5.17p and (|5.18p . Integrating (j5.12[) with respect to s and t, we obtain 

the second inequality in (|5.13|) . The first inequality in (|5.13p follows from the definition of 

and the monotonicity of r i— ► ^(pi >T ) (see the energy identity (|4.6p ). □ 

6 Proof of the main theorems 
6.1 The generation result 

Recall that §t{po) = St(po)^f d ^ when po = poJt? d \ n ; Theorem l4.11l is an immediate consequence 
of the following result. 
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Theorem 6.1. Let VL be a bounded, convex open set ofM. d and let us assume that /ij € M + (f2), 
i = 0,1, have finite distance W m> Q(/xo, fJ>i) < +00 (and therefore the same mass m = po(£l) = 
pi(Q)) and satisfy ! %f(po) < +°o and p\ = piJzf d | n <C ££ d ■ If U satisfies the generalized 
McCann condition GMC(m,d) (\4- TTp , h) then 

^ >n (S tA ti, tio) + tW(S m ) < ^W 2 |n ( Ml , fM Q ) + tW(fio) Vt > 0. 
Proof. Since W mj n(^o> /ii) < +00 there exists a geodesic curve (/x, 1/) 6 Q£q(pq — ► /zi) such that 

$m,n(Ms,»'s) = / ^m,n(Ai s ,^ s )ds = n (/X ,/Xl). 
JO 

Applying Lemma 13.51 we find a family of approximating curves {p e,5 ,v e,s ) and smooth convex 
open sets S\ satisfying (|3.7|) . (|3.8|) and 

JO 

Let m^P,, as in Lemma 021 for constants < M' < M" < M such that M' < p 6 ' 5 < M" in 
and let Cf' 5 '* 7 S C°°(M d ) obtained by solving (|5.3p with respect to and d s p £ > s in £ . Since 
m?7 > rn, by Theorem 5.21 of [DNS09J we easily have 

f \VC £ /' v \ 2 m v (p £ ' 5 ) dx = <S> mv ,nM' 5 ^ £ ' 5 ) < $m,n £ (/^ ,< V' 5 ) Vr/ > 0, s E [0, 1]. 

If g 5 ' 7 ? = S(^, m^, f2 e ) is the semigroup associated with S(P V , Q £ ) and the corresponding integral 
functional ()5.13|) gives 

V m ^ £ (SfVi'Vo' 5 ) +t%(Sfyf ) < i /' / jVC^| 2 m r? (^)dxd S + t^(^) 

<f^n(Mi,/*o) + *^Wj'^ 

Passing to the limit as 77 [ (notice that the functions S^ ,v p E,s take their values in [M', M"]) we 
get 

where S e = §(^',m, S\) is associated with S(P,Q e ). We can then pass to the limit as <5 j 0: 
since — > p e in L°°(f2 e ) we immediately have 



-.2 

2 ' ' m,ii£\~tri) ru/ ' " ™ vtri; — ^ 

Finally as e J, we conclude, recalling Proposition 14.31 □ 
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6.2 Geodesic convexity 

The proof of Theorem 14.131 follows immediately from the generation result Theorem 16.11 and 
Theorem 14.101 if every measure fj, s of the geodesic curve is absolutely continuous with respect to 
«Sf d | fi (see also Remark I4.15p . On the other hand, this property is not known a priori, so we 
need a more refined argument. 

Proof. As in Proposition ^. 41 we set /i e s := n s *k £ , u\ := v s *k e and we denote by S £ = S(^, m, Q e ). 
By (|3.5p and the contraction property given by Theorem 14.111 in Q e we have 

W m ^(8 e t fi £ Si ,§ £ t fi £ S2 ) < W m ,n £ (/4,/4) < w m,nK,fe) = \si - s 2 |W m) nOo,/^i) (6.1) 
and by <pU)j) . Theorem 0TTU] and P~T5j) we have 

^(§1/4) < (1 - s)W(jj® + sW(»l) + | S (1 - a), (6.2) 
W^^Sf/^) < *CW) ~ inf ^) < *(^(A*i) " mf (6.3) 



where the second inequality in (16. 3D follows from Jensen's inequality. 

We choose now a countable set ^ dense in [0, 1] and containing and 1, a vanishing sequence 



(£fc)fceN and another vanishing sequence (Ek)k£N so that limfc| +00 t~^ l 5l k = 0. By compactness 
and a standard diagonal argument, up to extracting a further subsequence, we can find limit 
points p, s for s £ such that 

§tfcMs fc — 51 fis weakly as k f +oo. 
By (I6.ip and Proposition 13.21 we get 

W m ,n(/i S i, As 2 ) < |*i - s2|W mj n(Mo,Mi) Vsi^e^f. (6.4) 



(I6.3|) yields /io = MO; fix = Mi so that we can extend /} to a continuous curve (still denoted by 
/i) connecting /io and fj,% still satisfying (|6.4|) for every s\,S2 € [0, 1]. The triangular inequality 
shows that (16. 4p is in fact an equality and the curve fl is a constant speed minimizing geodesic. 
On the other hand, the lower semicontinuity of ^ with respect to weak convergence and (|6.2p 
yields 

< (l-*)«r(Ato) + s^(Mi) Vse^. (6.5) 



A further lower semicontinuity and density argument shows that (16. 5|) holds for every s € 
[0,1]. □ 
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Final remarks and open problems 



This paper is a first step towards the investigation of the geometry of spaces of measures metrized 
by W m> Q, the induced convexity notions for integral functionals and the corresponding generation 
of gradient flows with applications to various nonlinear evolutionary PDE's. 

Since a sufficiently general theory is far from being developed and understood, it is in some 
sense surprising that one can reproduce in this setting the celebrated McCann convexity result. 
On the other hand, many interesting and basic problems remain open: here is just a provisional 
list. 

- At the level of the distance W mj f2 only partial results on some basic properties (such as 
density of regular measures or necessary and sufficient conditions ensuring the finiteness 
of the distance), are known and a complete and accurate picture is still missing. 

- The situation is even less clear in the case of unbounded domains: in this paper we 
restricted our attention to the bounded case only. 

- The study of other integral functionals is completely open, as well as applications to 
different types of evolution equations, like scalar conservation laws or nonlinear fourth 
order equation of thin-film type, see the introduction of [D NS09] . 

- It would be interesting to study other metric quantities (e.g. the metric slope) and the 
pseudo-Riemannian structure (tangent space, Alexandrov curvature, etc) connected with 
the distance and the energy functionals, see |AGS05|. ICMV06j . 

- The regularization properties and asymptotic behaviour of the gradient flow ^ and its 
perturbation can be studied as well: in the Wasserstein case the geodesic convexity of a 
functional yields many interesting estimates. 

- The convergence of the so called "Minimizing movement" or JKO-scheme could be ex- 
ploited in this and other situations: in the case of geodesically convex energies, fur- 
ther information on the Alexandrov curvature of the distance W m n would be crucial, 
see lAGS05ll5a07j . 
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